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ABSTRACT 

RX J0720.4-3125 is an isolated neutron star that, uniquely in its class, has shown changes in its thermal 
X-ray spectrum. We use new spectra taken with Chandra's Low Energy Transmission Grating Spectrometer, 
as well as archival observations, to try to understand the timescale and nature of these changes. We construct 
lightcurves, which show both small, slow variations on a timescale of years, and a larger event that occurred 
more quickly, within half a year. From timing, we find evidence for a 'glitch' coincident with this larger event, 
with a fractional increase in spin frequency of ~5 x 10~ 8 . We compare the 'before' and 'after' spectra with 
those from RX J 1308. 6+2 127, an isolated neutron star with similar temperature and magnetic field strength, but 
with a much stronger absorption feature in its spectrum. We find that the 'after' spectrum can be represented 
remarkably well by the superposition of the 'before' spectrum, scaled by two thirds, and the spectrum of 
RX J 1308. 6+2 127, thus suggesting that the event affected approximately one third of the surface. We speculate 
the event reflects a change in surface composition caused by, e.g., an accretion episode. 
Subject headings: stars: individual: (RX J0720. 4-3125, RX J1308. 6+2127) — stars: neutron — X-rays: stars 
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1. INTRODUCTION 

Among nearby neutron stars, seven show predominantly 
thermal emission, with inferred temperatures of kT = 50 
to 100 eV. These so-called isolated neutron stars (INS) are 
thought to be young, < 1 Myr old, with the emission due to 
residual heat. Compared to similarly aged pulsars, they have 
much weaker non-thermal emission, longer, 3 -10 s periods , 
and stronger, few 10 13 G mag netic fields (see HaberU 120061: 
Van Kerkwijk & Kaplan 2006 for recent reviews). 

The INS have attracted much attention, in part because of 
the hope that the equation of state in the ultra-dense interior 
could be constrained using properties derived from their spec- 
tra, such as effective temperature, surface gravity, and grav- 
itational redshift. Progress has been stymied, however, by 
difficulties in interpreting the spectra: at present, hydrogen, 
helium, and mid-Z elements, in states ranging from gaseous 
to condensed, all being cons idered (see, e.g., contributions to 
iPage. Turolla. & Zanell2006l) . 

A source that may elucidate matters is RX J0720.4-3125. 
This so urce was discover ed with ROSAT, and has a 8.39 s 
period (lHaberl et al.ll 1997b . Initial XMM-Newton and Chan- 
dra spectra appeared to be featurel ess, blackbody-like 
dPaerels etal.1 12001: Kap lan et all 12003), simil ar to those of 
the brightest INS, RX J1856.5-3754 (e.g., iBurwitz et ail 
12001b . but more detailed analys is uncovered a weak , broad 
absorption feature near 280 eV (HaberletalJ [2004a). sim- 
ilarly broad (but stronger) features at e nergies from 200 to 
700 eV were found in five other INS (lHaberl et all 120031 
I2004bt IVan Kerkwijk et alJ l2004t IZane et al l 120051) and are 
generally interpreted as proton cyclotron or neutral-hydrogen 
transitions. For RX J0720.4-3125, either interpretation is 
consistent with the magnetic field strength B~2x 10 13 G in- 
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ferred f rom timing (Kapl an & van Kerkwiikl2005al hereafter 
lKvK05t) . 

Unlike other INS, ho wever, the spectr um of 
RX J0720.4-3125 changed dDe Vries et al.1 l200l . The 
change could be described as an increase in temperature, 
from 86 to 94 eV, and an increase in equivale nt width of 
the a bsorption feature, from ~ 5 to 70 eV (Hab erl et al.l 
2006, hereafter lH06l) . The nature of the cha nges is unclear. 
De Vries et ail d2004l) and IVink et al.l (120041) find that they 



started gradually and then speeded up, and lH06l find a recent, 
partial reversal. These authors interpret this as evidence for 
precession on a ^7yr period. 

In this Letter, we reconsider the timescale of the changes, 
investigate what changed, and offer an alternative interpre- 
tation. After describing the data we use in §|2] we present 
lightcurves in §[3] and update the timing solution of iKvKOBl 
in §H] finding evidence for a glitch at the time of the largest 
spectral change. We study the 'before' and 'after' spectra in 
§[5] and suggest that only part of the surface changed. We 
discuss implications in §|S] 

2. OBSERVATIONS, REDUCTION, AND ANALYSIS 

We used Chandra to observe RX J0720.4-3125 several 
times in 2005 and 2006, for ^300ks total, taking spectra with 
the Low Energy Transmission Grating Spectrometer (LETG). 
For comparison, we analysed archival spectra taken in 2000 
(~38ks) and 2004 (^36ks), as well as spectra of another 
INS, RX J1308.6+2127, taken in 2004 (~90ks). 

For all observations, we extracted spectra following stan- 
dard procedures (we used CIAO ver. 3.3 and CALDB ver. 
3.2). We calibrated the spectra outside of CIAO: We sub- 
tracted the background, corrected for higher-order contribu- 
tions using the observed counts and the Nov. 2004 effective 
areas, 4 and calculated photon rates using the first-order effec- 
tive areas. Note that we implicitly assume diagonal response 
matrices, thus ignoring the LETG line-spread function; this 
should be fine for spectra without sharp features, like ours. Fi- 
nally, for timing purposes , we extracted barycentered source 
events (following |KvK05l) . 

4 |http : / / cxc . harvard ■ edu/cal/Letg/Hrc_QE/ea_index .htmT| 
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FIG. 1 . — Lightcurves and timing residuals for RX J0720.4— 3125. Panels 
a and b: Photon rates in various wavelength bands. The central wavelengths 
are indicated; the widths are 5 A for bands between 15 and 35 A, and 10 A 
for those from 45 to 65 A (the latter bands are not covered by RGS). For clar- 
ity, the 30 A fluxes are offset vertically by 0. 1 unit. The formal uncertainties 
are smaller than the symbol size. Panel c: Arrival-time residuals relative to 
a simple spin-down model (fit to Chandra and XMM data only). The differ- 
ent instruments are indicated, and the uncertainties include, in quadrature, a 
systematic term of 0.1 1 s. Overdrawn (dashed red curve) is the best-fit glitch 
model, which reproduces the data much better. Panel d: Difference in spin 
frequency of the glitch model relative to the simple spin-down model (solid 
curve with lcr uncertainties indicated in grey). 

We augmented our data with archival XMM observations. 
For timing purposes, we extracted ba rycenter ed source events 
from all EPIC instruments (following KvK05), while for mak- 
ing lightcurves we extracted RGS spectra. The latter was done 
with rgsproc (XMMSAS ver. 7.0), which was also used 
to produce flux-calibrated spectra (ignoring, like for LETG, 
the line-spread function). We note that XMMSAS ver. 7.0 
includes corrections for the changing long-wavelength sensi- 
tivity of RGS, 5 which has influenced previous analyses (see 
§0. 

3. TIME SCALES OF THE SPECTRAL CHANGE 

To determine the timescale on which the spectral changes 
occurred, we used the LETG and RGS spectra to produce 
lightcurves in eight wavelength bands (Fig.[TJ. Broadly, these 
confirm previous results (see §[TJ: as time passed, the flux in- 
creased at wavelengths A < 25 A, and decreased in the range 
35 < A < 55 A, presumably because the absorption feature 
deepened. However, with the new RGS long-wavelength ef- 
ficiencies, the evidence for a decrease in 25-36 A flux before 
2003 has largely disappeared (but the first LETG and RGS 
observations have become discrepant, suggesting some cali- 
bration issues remain). 

Figure[T]shows that most of the change occurred rapidly, in 
between the 2003 May and 2003 October XMM observations. 

5 See XMM-SAS ver. 7.0 and XMM-CCF-REL-216 release notes. 



TABLE 1 

Log of New Observations and Times of Arrival 

Exp. TOA c 
Instrument" ID b Date (ks) Counts (MID) 

PN 986-S3 2005 Apr 28 52.0 300204 53488.6756038(3) 

MOS1 986-S1 41.1 51871 53488.6013068(8) 

MOS2 986-S2 41.2 54046 53488.6012096(8) 

HRC 5582 2005 lun 01 69.7 903822 53522.9397230(12) 

6364 2005 Aug 27 38.6 520744 53610.0880958(18) 

PN 1060-S3 2005 Sep 23 51.0 284343 53636.3001458(4) 

MOS1 1060-S1 36.9 50393 53636.2024428(8) 

MOS2 1060-S2 36.9 51231 53636.2027347(8) 

HRC 6369 2005 Oct 08 26.0 409201 53652.0859352(21) 

7177 2005 Oct 09 8.0 106581 53652.5397755(32) 

PN 1086-S3 2005 Nov 12 38.0 268168 53687.1718885(4) 

MOS1 1086-S1 39.5 56034 53687.1658660(8) 

MOS2 1086-S2 39.5 56292 53687.1658654(7) 

HRC 7243 2005 Dec 14 17.1 241644 53718.1462899(25) 

7244 2005 Dec 15 16.2 227541 53719.8089713(22) 

7245 2005 Dec 16 17.1 249521 53720.7701604(30) 
5584 2005 Dec 17 14.1 194582 53721.9195718(30) 

7251 2006 Feb 09 10.6 160780 53775.3508981(32) 

NOTE. — A log of earlier observations and corresponding times of arrival 
can be found in Kaplan & van Kerkwiik (2005a). 

a HRC: High-Resolution Camera for spectroscopy aboard Chandra, used with 
the Low-Energy Transmission Grating. PN: European Photon Imaging Cam- 
era (EPIC) with PN detectors aboard XMM, used in full-frame mode with a 
thin filter. MOS1/2: EPIC with MOS detectors aboard XMM, used in small- 
window mode with thin filter. b Observation identifier (Chandra) or revolution 
number and exposure identifier (XMM). C Time of arrival of maximum light 
closest to the middle of an observation, as determined from fitting a sinusoid 
to the binned lightcurve; it is given with l-cr uncertainties. 

Before and after, there are slower variations, with those before 
in the same sense as the rapid change, and those after in the 
op posit e one. This is consistent with inferences from PN data 
bvlH06l 

To examine broad-band changes, we integrated the LETG 
fluxes in the 10-80A range (0.155-1. 25 keV), finding totals 
of 0.97, 1.06, 1.04, 1.02, 1.04, and 1.03 x 10" 11 ergs" 1 cm" 2 
for the six observations, respectively (with — 2% uncertainty). 
Thus, from 2000 to 2004 the flux increased by -9%, and af- 
terwards it possibly decreased slightly. The change of the total 
source flux is likely similar: most flux is emitted in the 10- 
80 A band, and, despite the change in spectral shape, different 
epochs have extinction corrections identical to < 1% (e.g., 
factor 1 .43 for N H = 1 .3 x 10 20 cnT 2 ). 

H06 suggested the changes were cyclic, with a ~ 7 yr pe- 
riod, but the lightcurves do not require this. Rather, to us 
Fig.[T](as well as Fig. 2 of H06) suggests a major, impulsive 
change in 2003, the source 'building up' to the change before- 
hand, and 'relaxing back' afterwards. 

4. A CHANGE IN TIMING PROPERTIES 

We updated the timing solution of lKvK05l determining ar- 
rival times for all new observations (see TableQ}. These were 
consistent with the old timing solution, but, as before, signif- 
icant timing residuals remained (Fig. [TJ: to obtain x 2 ed = 1 
required adding (in quadrature) a systematic error term of 
-0.3 s. 

Inspired by what to us seems a sudden spectral change (§0, 
we considered a 'glitch' as the source of the timing residuals. 
Earlier, we had found some improvement with a rudimentary 
glitch model, in which the freque ncy v c hanged after mid- 
2003 and recovered only slightly (lKvK05h . Here, we tried a 
different model, in which we fit for different v and v before 
and after a given glitch time t g , but keep phase continuous, ig- 
nore any recovery, and leave the cycle counts unchanged from 
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TABLE 2 

Timing Solutions for RX J0720. 4-3125 



(Quantity 




All Pioto 


Spindown only 






to (MJD) 


53010.2635646(7) 


53010.2635626(6) 


v (Hz) 


0.11917366979(12) 


0.11917366954(11) 


!/(Hzs _1 ) 


-9.74(4) x lO" 16 


-9.88(13) x 10~ 16 


TOA rms (s) 


0.26 


0.29 


X 2 /dof 


77.6/46=1.69 


150.8/49=3.08 


Spin-down + Glitch 






t (MJD) 


53010.2635686(10) 


53010.2635667(10) 


y(Hz) 


0.1191736716(9) 


0.1191736716(9) 


i>(10~ 15 Hzs-') 


-1.04(3) 


-1.04(3) 


t g (MJD) 


52817(61) 


52866(73) 


At/ (nHz) 


5.7(17) 


4.1(12) 


Ai>(l(T 17 Hz s _I ) 


-1(4) 


-4(3) 


TOA rms (s) 


0.15 


0.24 


X 2 /dof 


37.0/43=0.86 


45.1/46 = 0.98 



NOTE. — The parameters determine the cycle count plus phase 
via 0(0 = v(t -t )+ ^v(t -t ) 2 + A<f)g(t), where Arf> g (t) = -Au(t - 
t g ) — j Ai>(t — t g ) 2 for t < t g in the glitch model and zero other- 
wise. For all fits, a 0. 1 1 s systematic uncertainty has been added 
in quadrature to the times of arrival (TOAs), and the uncertainties 
quoted are twice the formal lcr values. 

the simple spin-down model. This improves the fit drastically 
(Fig-tD: to get x 2 ed ~ 1, a systematic error term of only 0.1 1 s 
is required. Intriguingly, the best-fit glitch time coincides with 
the period of the greatest spectral change. 

In Table|2] we list parameters for the simple spin-down and 
the glitch models. We try solutions both with and without the 
ROSAT data (which have much less secure cycle counts); in 
either case, the change in v is significant, but that in v is not. 

The addition of a glitch improves the fit, but also adds com- 
plexity. To evaluate the significance of the improvement, we 
tried two other ways of adding complexity. First, we in- 
cluded higher-order frequency derivatives, but found a sim- 
ilarly improved fit required eighth order (i.e., six ad dition al 
parameters). Second, motivated by the suggestion of lH06l of 
a ~ 7 yr periodicity in the spectral variations, we added a pe- 
riodic component, but found a similarly improved fit only for 
periods near 4yr (as expected from Fig. [T] for details, see 
|VanKerkwijk & Kaplan 2006). The result did not seem to 
depend on whether or not we included the ROSAT data o r 
changed their cycle counts (as suggested by lHaberll 120061) . 
Overall, the timing seems to support an impulsive model, but 
does not exclude periodic changes. 

5. NATURE OF THE SPECTRAL CHANGE 

The lightcurves show that most of the spectral change in 
RX 10720.4-3125 took place in 2003. In Fig. El we show 
the LETG spectrum from 2000, taken before this time, as 
well as the average of all LETG spectra taken after 2003 (the 
further small changes in 2004 and 2005 [Fig. [T| would not 
be noticeable in Fig. |2). For comparison, we also show the 
LETG spectrum of RX 11308.6+2127, which is very similar 
to RX 10720.4-3125 in its inferred temperature, extinction, 
and magnetic field st rength (kT ~ 1 eV, jV H - 2 x 1 20 cm" 2 
Schwope et ail20 06: B-3x 10 13 G. lKaplan & van Kerkwijkl 
20059) - 




The comparison of the spectra raises two main questions. 
First, for RX 10720.4-3 125, how can the appearance of a pro- 
nounced absorption feature be associated with o nly a small 
change in temperature (from kT ~ 86 to 94 eV; H06)? Sec- 
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FIG. 2.— Calibrated LETG spectra. The spectrum of RX J0720.4-3 125 
from 2000 (red open circles) is well described by an absorbed black-body 
(jV H = (1.55 ±0.11) x 10 20 crrT 2 , kT = 83.3 ± 1.1 eV, and r/d = 0.0189 ± 
0.0008 km pc~' ; drawn red curve). The average post-2003 spectrum (ma- 
genta filled circles) could be a superposition of spectra from unchanged and 
changed parts of the surface. As an indication, a close to maximum contri- 
bution by the unchanged spectrum is indicated (magenta, dotted line). The 
residual (blue, wiggly curve, offset down by 0. 1 unit for clarity) is remark- 
ably similar to the spectrum of RX J1308. 6+2127 (green squares, scaled by 
a factor 1.2). For all spectra, the 40-44 A range is less reliable, because of 
strong instrumental carbon absorption. 

ond, why is the absorption feature in RX 11308.6+2127 much 
stronger than that in RX 10720.4-3125 (in either state)? 

Of these two questions, we believe the answer to the sec- 
ond may be simple: that not the whole but only part of the 
surface changed. If so, then the 'after' spectrum contains a 
contribution from unchanged parts that emit the cooler, nearly 
featureless 'before' spectrum. Hence, the spectrum from the 
changed part should be hotter and should have a stronger line 
than one would infer from the average. 

We do not know the relative contribution of the 'before' 
spectrum, but an upper limit of ~ 70% is set by the re- 
quirement that it does not exceed the 'after' spectrum at 
any wavelength (set by the 35^-0 A region; Fig. |2). For 
a contribution near this limit, the absorption feature in the 
changed parts of the surface would be very strong, just as ob- 
served in RX 11308.6+2127. Indeed, if we subtract the 'be- 
fore' spectrum, scaled by two thirds, from the 'after' spec- 
trum, the residual is remarkably similar to the spectrum of 
RX 11308.6+2127 (see Fig. |2). This is confirmed by the 
EPIC-pn data (Mori et al., in preparation). 

6. RAMIFICATIONS 

We have found that most of the spectral variations in 
RX 10720.4-3125 occurred on a short, less than half-a-year 
timescale, and that, coincident with this, there likely was an 
increase in spin frequency. We also showed that the 'after' 
spectrum could be reproduced if one third of the emission 
changed from the blackbody-like 'before' spectrum to some- 
thing similar to the emission of RX 11308.6+2127, and we 
suggested only part of the visible surface changed. Here, we 
will interpret these findings assuming an impulsive event (see 
|H06| for a discussion in terms of precession). As the obser- 
vations show evidence for a reversal, we need to address not 
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only the nature and cause of the change, but also how it can 
be undone. 

It is difficult to understand the change without knowing the 
composition and state at the photosphere (§[TJ- The order-of- 
magnitude increase in strength of the absorption feature sug- 
gests, however, that a new species became dominant, due to 
a change in ionisation balance, state (e.g., from condensed to 
gaseous), or composition. 

In principle, a temperature change could affect the ionisa- 
tion balance or state, and one could envisage heating due to a 
sudden energy release associated with, e.g., a glitch, and sub- 
sequent cooling on a timescale set by the release depth. How- 
ever, the average observed flux before and after the event dif- 
fers by only ~5%. Even if due to only one third of the surface, 
the im plied increase in effective temperature is at most 4% 
(H06 find a 6% change in average temperature from spectral 
fits, but this may simply reflect that a model of a blackbody 
with a gaussian absorption line is not realistic.) With such 
a small change, it is not clear the atmosphere and spectrum 
could change so drastically. 

An easier explanation for the strengthening of the absorp- 
tion feature would be a change in composition, e.g., due to 
an accretion episode. Such a change could also lead to a 
small increase in effective temperatur e, since heat conduc - 
tion is faster for lighter elements (e.g., Pote khin et al.ll2003h . 
But how to undo the change? For hydrogen, ther e is a pos- 
sible mechanism (Chang . Arras. & Bildstenl [20041) : for B ~ 
2 x 10 13 G and T eff ~ 10 6 K, hydrogen will be burnt diffu- 
sively on carbon in ^200y H 5 (o 2 yr (where = 10)?H,iogcm~ 2 
is the hydrogen column density), with substantial uncertainty 
due to the composition of the underlying layers. A thin hy- 
drogen layer , superposed on a condensed su rface, has been 
invoke d by iMotch. Zavlin. & Haberll (120031) and iHo et al l 
(2006) to explain the featureless spectra of some INS. For 
RX J0720.4-3125, it would imply a very thin hydrogen layer 
before the change, and a thicker one on parts of the surface 
afterwards. 

Among the two possibilities mentioned above - a glitch 
in the neutron-star crust and an accretion episode - the for- 
mer would most naturally explain the timing results; in- 
deed, the fractional increase in spin frequency of Av/v = 
5 x 10~ 8 would not be unu sual for an older radio pulsar (e.g., 
iJanssen & Stapperi f2006). Also for magnetars, glitches are 



observed, as well as slow changes on timescales of years (see 
lKaspill2006l for a review). The glitch could release some en- 
ergy, causing a t emperature increase that depends on the re- 
lease depth (e.g.. lHui & Chenal2004l and references therein). 
In addition, it could induce magnetospheric currents, which 
could produce hydrogen and helium on the surface by spal la- 
tion (as suggested for RX J 1856.5-3754 by |Ho_et_aL||2006). 

Accretion would naturally lead to a change in composition, 
as well as, by its associated angular momentum, a change in 
spin. Matter could originate in an asteroid or comet belt such 
as has been invoked t o explain timing no i se, n ulling and other 
pulsar properties by ICordes & Shannon! (120061) . If accretion 
started due to a collision between asteroids, low-rate accre- 
tion of smaller debris could perhaps cause the slower changes. 
The observed frequency change requires an accreted mass of 
~ 10 20 — 10 21 g, depending on accretion regime and properties 
of accreting matter (large bodies accreting directly onto the 
star, small grains destructed and ionized before reaching the 
Alfven surface, etc.). The accretion would result in an en- 
ergy of ~ 10 40 -10 4I erg being released in a relatively short 
time, hence a high luminosity. We have checked the RXTE 
All-Sky Monitor (ASM) data in 2003 May-October and found 
no emission from RX J0720.4-3125; however, given the low 
ASM duty cycle, a short-lived (< lOks) event could easily be 
missed. 

In future, the above ideas - a s well as the alterna- 
tive suggestion of precession (lH06l) - can be tested. For 
instance, both give different predictions for future tim- 
ing. Generally, if hydrogen appeared on the surface (or 
hydrogen-covered parts became visible), the optical flux 
should perhaps have increased, yet it ch anged by < 10% 
(Kapl anTvan Kerkwiik. & Andersonll2007l) . Also, if the com- 
position changed locally, e.g., near a magnetic pole, mat- 
ter might diffuse outwards, perhaps leading to further slow 
changes. If so, one would also expect that the pulse pro- 
file evolved, perhaps being most sharply peaked, especially 
at high energies, shortly after the event. From the PN data, 
this indeed appears to be the case (Mori et al., in preparation). 

We thank Ron Remillard for help with the RXTE ASM 
data analysis. Support for this work was provided by NASA 
(Chandra award GO5-6074). 
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